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Abstract 

We consider radiative electromagnetic corrections, at order a, to the process e + e _ —> 
vV at finite density and temperature. This process represents one of the main con¬ 
tributions to the cooling of stellar environments in the late stages of star evolution. 
We find that these corrections affect the energy loss rate by a factor (—4 4 - 1)% 
with respect to the tree level estimate, in the temperature and density ranges where 
the neutrino pair production via e + e _ annihilation is the most efficient cooling 
mechanism. 
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1 Introduction 


The late stages of star evolution are strongly influenced by neutrino emis¬ 
sion processes. After hydrogen burning, stars with small mass, of the order 
of the solar mass, evolves towards a white dwarf configuration. On the other 
hand more massive stars undergo several burning phases till the formation of 
a central iron core and the eventual triggering of the gravitational instabil¬ 
ity leading to the Supernova phenomenon. In both cases the cooling rate is 
largely dominated by neutrino production. For example, from the 12 C — 24 Mg 
burning phase, almost 100% of the energy flux is emitted via neutrino produc¬ 
tion, which leave the stellar system without any interaction because of their 
extremely large mean free path. A precise determination of neutrino emission 
rates is therefore a crucial issue in any careful study of the final branches of 
star evolutionary tracks. In particular, changing the cooling rates at the very 
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last stages of massive star evolution may sensibly affect the evolutionary time 
scale and the iron core configuration at the onset of the Supernova explosion, 
whose triggering mechanism is still lacking a full theoretical understanding [1], 

The energy loss rate due to neutrino emission, hereafter denoted by Q, receives 
contribution from both weak nuclear reactions and purely leptonic processes. 
However for the rather large values of density and temperature which char¬ 
acterize the final stages of stellar evolution, the latter are largely dominant. 
They are mainly due to four possible interaction mechanisms: 


i) 

pair annihilation 

e + + e - 

-> v + V 

ii) 

z/-photoproduction 

7 + e ± - 

-> + v + 77 

hi) 

plasmon decay 

* 

7 - 

-> V + V 

iv) 

bremsstrahlung on nuclei 

e ± + Z - 

-> e ± + Z + v + 77 


Actually each of these processes results to be the dominant contribution to Q 
in different regions in a density-temperature plane. For very large core tem¬ 
peratures, T > 10 9 °K, and not too high values of density, pair annihilations 
are most efficient, while v photoproduction gives the leading contribution for 
10 8 °K < T < 10 9o /l and relatively low density, p < 10 5 g cm -3 . These 

are the typical ranges for very massive stars in their late evolution. Finally, 
plasmon decay and bremsstrahlung on nuclei are mostly important for large 
(p ^ 10 6 g cm -3 ) and extremely large (p > 10 9 g cm -3 ) core densities, re¬ 
spectively, and temperatures of the order of 10 8 °K < T < 10 10 °K. Such 
conditions are typically realized in white dwarfs. 

Starting from the first calculation of neutrino energy loss rates, based on the 
V-A theory of weak interactions [2,3], a systematic study of processes i)-iv) 
has been carried out in a long series of papers [2]-[22]. In all these analyses 
the pair production rate i) has been evaluated at order Gp, i.e. at the zero 
order in the electromagnetic coupling constant a expansion. Despite of the 
different topology and phase space volume which characterize the remaining 
processes ii)-iv), it is nevertheless worth noticing that they are instead at least 
of order aG 2 F . In view of this, it is therefore meaningful to investigate whether 
including QED radiative corrections to pair annihilation rate i) may lead to 
a sensible change in the cooling rate Q. This is the aim of this letter. In par¬ 
ticular we report here the results of a calculation at order a of two classes 
of contributions, due to vacuum radiative corrections as well as those arising 
from the interaction of the e + -e“ pairs with the surrounding electromagnetic 
plasma. This has been performed in the real time formalism for finite temper¬ 
ature held theory. We actually fold that the total corrections to Q range in the 
interval (—4 -7 1)% of the tree level estimate, in the temperature and density 
ranges where the pair annihilation process represents the main contribution 
to the total energy loss rate. We have also re-evaluated both photoproduction 
and plasmon decay rates, and we do find a good agreement with the most 
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recent estimates reported in the literature [21,22], 

The paper is organized as follows. Section 2 is a short review of the lowest 
order calculation of Q for the pair annihilation process. The order a QED 
corrections are then considered in Section 3, where we report the main results 
of our calculations of both vacuum and thermal contributions. We finally 
discuss our results and give our conclusions in Section 4. 


2 Pair annihilation in Born approximation 


We first consider the pair annihilation process e + + e~ —> v + v at the lowest 
order in perturbation theory. Denoting with pig = (Tag,Pi, 2 ), the energy- 
momentum of ingoing electron and positron, and with gig = (uqg, gig) the 
four-momentum of the outgoing neutrino and antineutrino, respectively, the 
neutrino energy loss rate Q is given byQ 


n 1 r d 3 pi r d 3 p 2 

J ~ (2tt) 6 J 2E 1 J 2 E 2 


(E 1 + E 2 )F_(E 1 )F + (E 2 ) 



( 1 ) 


where 


F±(E) = exp {f ±£ e } + 1 


-1 


are the Fermi-Dirac distribution func¬ 


tions for positrons/electrons with temperature T and degeneracy parameter 


Finally £ a ,/3 |M e+e -_ 


i2 • 


I 


is the squared modulus of the reaction ampli¬ 


tude summed over all particle polarizations A and neutrino flavours (5. 


In the Born approximation, i.e. in the limit of a four fermion electroweak inter¬ 
action and no electromagnetic radiative corrections (see Fig. la), the quantity 
in curly brackets in Eq. (1) takes the form 


2 G 2 f 

37T 


{pv + C'a ) (m 2 e + 2p r p 2 ) + 3 (Cy 


cr 


) P] (P1+P2Y 


( 2 ) 


where we have defined Cy A = (1 + Cv,a ) 2 + 2 Cy A with Cy = 2 sin 2 9\y — 1/2, 

C A = —1/2. 

Substituting Eq. (2) into Eq. (1) and performing the angular integrations, we 
get the final expression for the energy loss rate induced by pair annihilation 
in the Born approximation, 

1 We use natural units, h = c = k = 1. 


3 
















Fig. 1. Feynman diagrams for the pair annihilation process up to order aG 2 F . 
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E 2 + E 2 | { E 2 E 2 \ 
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( 3 ) 


Note that Qb depends on the temperature T and the electron degeneracy 
parameter £ e only. It is customary to express the dependence on £ e (or the 
electron chemical potential) in terms of the matter density, p, the electron 
molecular weight, p e , and T through the condition of plasma electrical neu¬ 
trality, 


P_ 

pe 



N J 4 being the Avogadro number. 

In Fig. 2 we show the energy loss rate due to pair annihilation in the Born 
approximation for several values of the temperature as function of p/p e . 
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Fig. 2. The energy loss rate due to pair annihilation process in the Born approxi¬ 
mation for T = 10 s , 10 8-5 ,10 9 ,10 10 °K (from bottom to top). 

3 Pair annihilation: radiative corrections 


A consistent computation np to order a of the energy loss rate Q due to 
pair annihilation is obtained by considering the additional contribution of all 
diagrams shown in Fig.s lb-lh. Notice that since the energy carried by the 
outgoing neutrino pair is typically smaller or, at most, of the order of 1 MeV, 
we can safely neglect the electroweak radiative corrections to the four fermion 
effective interactions, and only consider the set of gauge invariant purely QED 
contributions. 

The radiative corrections in case of large values of temperature and density 
must be computed taking into account the presence of the electromagnetic 
plasma, namely e ± and 7 . As long as neutrino mean free path is so large 
that they leave the star without any further interaction, there is no relevant 
neutrino component in the stellar plasma. The zero temperature radiative 
corrections, or vacuum radiative corrections, have been already computed in 
literature. They can be obtained for example, from the results of Ref. [23] by 
using crossing symmetry. It is customary to classify them as follows: 

- electron mass and wavefunction renormalization (Fig.s lb-lc); 

- electromagnetic vertex correction (Fig. Id); 

- 7 emission/absorption (Fig.s le-lh). 

The whole radiative corrections, including the finite temperature effects, can 
be then obtained from the diagrams of Fig.s lb-lh by using the real time 
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formalism for finite temperature quantum field theory. In this framework the 
thermal propagators for electrons and photons read as follow 


St(p) = O + m){^— 


p* — rni + ie 


-2nd (p 2 -mfj [F_(|p o |) 0 (po) + i r +(bo|) 0 (-Po)]} 

Df(k) = - 


k 2 + ie 


+ 2nS(k 2 )B(\k 0 \) g<# 


( 5 ) 

( 6 ) 


where (~)(x) is the step function and B(x) is the Bose distribution function. The 
first terms in Eq. (5) and ( 6 ) are the usual T = 0 Feynman propagators, while 
those depending on the temperature (and density) through the distribution 
functions describe the interactions with real particles of the thermal bath. 


The calculation for the order aG 2 F corrections to the Born approximation for 
the energy loss rate proceeds in close analogy with what has been described in 
Ref.s [24,25] for a very similar situation. In particular, the additional contri¬ 
butions to Q due to mass, wave function and vertex corrections come from the 
interference between diagrams of Fig.s lb, lc and Id with the Born amplitude 
of Fig. la. In addition, the energy loss rates due to 7 emission (absorption) 
Q e (a), are given by the sum of squared amplitudes of the processes of Fig.s le, 

If (lg, lh)> 



For the seek of brevity we do not report here the details of our calculation leav¬ 
ing this lengthy description to a forthcoming paper [26]. We only briefly outline 
the method adopted to obtain the final corrections A Q. Each of the several 
contribution to \M e + e -^. VpVp | 2 , and the two squared amplitude \M e + e -_> VpVpl \ 2 
and \M e + e - 1 ^ VpVp \ 2 have been evaluated analytically. The results should be 
then integrated over the relevant phase space. Some of these integrations can 
be again done analytically, while we have used a Montecarlo technique to deal 
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Fig. 3. The total radiative corrections normalized to the Born approximation re¬ 
sult for the pair annihilation process for T = 10 s , 10 8 ' 5 ,10 9 ,10 10 °K (from top to 
bottom). 

with the (usually three of five-dimensional) remaining integrations, requiring 
an accuracy better than 1% on A Q 2. The main difficulty of this procedure 
is that each contribution coming from the interference of diagrams lb)-ld), 
as well as the squared amplitudes le)-lh) is plagued by infrared divergen¬ 
cies, while the overall result is of course divergence free. To regularize these 
divergencies we have explicitly subtracted all divergent terms expanding the 
squared amplitudes in a Laurent series around the pole singularities. This 
method follows quite closely what has been already used in [24,25,27]. 

In Fig. 3 we show the ratio between the sum of all radiative corrections and 
Qb , A Q/Qb, for some values of T as a function of the density. For a fixed 
value of temperature, the ratio A Q/Qb, shows a plateau, for all densities for 
which the electron-positron chemical potential is smaller than the tempera¬ 
ture, so that the mean energy is essentially given by T. The plateau end point, 
as expected, shifts towards larger values of p^ e with increasing temperature. 
For even larger densities the mean energy of the e + — e~ pair is dominated by 
their Fermi energy tp, and the ratio A Q/Qb start varying logarithmically as 
a function of tp/wie- Similarly, at fixed density and increasing temperature, 


2 The accuracy is actually depending on the values of temperature and density. The 
estimate of 1% for the theoretical error is rather conservative, and represents the 
larger result in the whole interval where pair annihilation significantly contribute 
to the cooling rate Q, 10 9 °K < T < fO 10 °K . For smaller values of T, since the 
pair annihilation process rate is very small, the accuracy gets worst, and in our 
calculations can reach 5%. 
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where the effect of degeneracy can be neglected, the ratio A Q/Qb shows a 
similar logarithmic behaviour as function of T/m t \ Both these behaviours 
suggest that in these extreme regimes higher order QED contributions should 
be included. However, the trends exhibited in Fig. 3 cannot be straightfor¬ 
wardly extrapolated where T or start to be closer to the W and Z masses, 
since the effects of electroweak corrections should be taken into account as 
well. Furthermore it is also worth noticing that for temperatures larger than 
10 10 °K, and/or densities larger than, say 10 12 g cm -3 , the neutrino mean free 
path becomes of the order or smaller than the dimension of the star core, so 
that neutrinos start interacting with the medium. This means that to eval¬ 
uate the cooling rate via neutrino emission it is necessary to consider their 
Pauli blocking factors, as well as the (negative) contribution due to inverse 
processes. 

It is actually interesting to see how the increasing importance of radiative 
corrections with temperature and density in the range of interest emerges 
from the study of other parameters which characterise the properties of the 
electromagnetic plasma, like thermal effective masses of particles. In general, 
thermal masses would show an explicit dependence on the spatial momentum 
p. Nevertheless a nice way to see how their dispersion relations are influenced 
by interactions is to consider their average over the equilibrium distribution. 
In particular, for the electron effective mass we have 

1 00 

\p\ 2 d\p\mfm)F_(E) , (9) 

7 t- Tl- J 
0 


where n_ is the electron number density and m^(|p|) is the momentum- 
dependent renormalized mass in the electromagnetic plasma with temperature 
T and density p//i e , 


d,, .. an T 2 a 

= ™ e + —-+- 

3 m P 7T m P 


\k\ 2 d\k\ . _ . . „ . .. am, 

11 11 (F4u) + F + (u)) 


\k\ d\k\ 


U) 


2n\p\ 


ui 


x F_ (a;) log 


ujE — ml + \p\\k\ 
ujE — ml — \p\\k\ 


- F + (u) log 


/ ujE + m 2 + 


\u)E + 


mt 


\m\ 


( 10 ) 


where a j — \]\k\ 2 + m 2 . An analogous result for positron can be obtained from 
(10) with the substitutions F±{u) —> Fn (a;). 


3 See, for example, the typical behaviour of QED radiative corrections for a similar 
process in [23] 
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Fig. 4. Contour plots for the photon (dashed lines) and the electron (solid lines) 
thermal mass correction, normalized to m e . The electron effective mass has been 
averaged over the thermal equilibrium distribution. From left to right the curves 
refer to the values 10~ 3 , 10 2 , and 10 _1 . 

For the photon thermal mass the averaging is unnecessary, since it does not 
depend on the spatial momentum 


m 


2 

7 


4a J \k\ 2 d\k\ 

IT J LU 
0 


(F.(u) + F + {u)) 


( 11 ) 


In Fig. 4 we have reported the contour plots of (< mf > — m e )/m e and 
m 7 /m e in the temperature-density plane. Curves from left to right refer to the 
values 10 -3 , 10 -2 and 10” 1 , respectively. From this plot it is clear that both 
electron and photon effective masses grows logarithmically with the particle 
mean energy, and are again of the order of few percent in the relevant range 
of temperature and density. 


4 Results and conclusions 


In this letter we have reported on the calculation of the energy loss rate of 
pair annihilation process i) up to order aG\. The Born value for Q has been 
corrected including both vacuum and thermal radiative corrections, the lat¬ 
ter being computed in the real time formalism. One of our main results is 
presented in Fig. 3, where we plot these radiative corrections with respect to 
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the Born approximation estimate, as functions of the plasma density for some 
values of the temperature. The corrections come out to be of the order of few 
percent and negative for high temperatures, implying that for these temper¬ 
atures the energy loss is sensibly decreased. At fixed temperature, A Q/Qb 
goes to a constant value for low density. This can be easily understood, since 
in this limit the plasma is weakly degenerate, and therefore the energy loss 
rate depends on temperature only. At large densities the ratio A Q/Qb de¬ 
creases and reach larger negative values. However for such high densities the 
pair annihilation rates are exceedingly small and thus this process gives only 
a marginal contribution to the star cooling. 

In order to perform a comprehensive study of order aG 2 F contributions to 
the energy loss rate, we have also recalculated the contribution from //-photo- 
production and plasmon decay processes. While the detailed analytical cal¬ 
culations will be discussed elsewhere [26], we here give our numerical results 
in Fig. 5, where we show the several contributions to Q. When possible, we 
also show for comparison the results of Ref. [22], In particular our results for 
these processes agree quite nicely, at least in the whole region where they sig¬ 
nificantly contribute to the total energy loss rate. For the bremsstrahlung on 
nuclei we have used the analytic fitting formula of Ref. [22]. 

Fig. 5 shows that, as well known, for large temperatures and not too high 
densities, pair annihilation dominates over the other two processes, while for 
low densities z/-photoproduction dominates over plasmon decay. On the other 
hand, for large densities the most relevant process is plasmon decay, whose 
rate however, along with those of all other processes, rapidly falls down for 
extremely high densities. This is a genuine plasma effect. Consider, for ex¬ 
ample, the behaviour of z/-photoproduction energy loss. As already noted in 
[3], the decrease for very large densities is achieved only if one consistently 
takes into account the increasingly large photon thermal mass. In fact with a 
massless photon the z/-photoproduction curves in Fig. 5 would rather reach a 
constant value. The main effect of m 2 is a lowering of the values of the Bose 
distribution function for photons, i.e. a smaller number of thermal photons. 
This reduces the energy loss rate induced by z/-photoproduction. 

In Fig. 6 we show the regions in the temperature-density plane where a given 
process contributes to the total energy loss rate (including radiative correc¬ 
tions to pair annihilation) for more than 90%. We also summarize there our 
results on the radiative corrections to pair annihilation processes, by plotting 
the contours corresponding to A Q/Q Fot = 1%, 0%, —1%, —2%, —3%, —4%, 
where Q Fot is the total emission rate with pair annihilation calculated in Born 
approximation. These contours lie almost entirely in the region where the pair 
annihilation process gives the main contribution to Q. This result may affect 
the late stages of evolution of very massive stars by changing their configura¬ 
tion at the onset of Supernova explosion. This issue is presently under study. 
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Fig. 5. The energy loss rate versus p/p e due to pair annihilation (including radiative 
corrections), photoproduction and plasmon decay (solid lines) for several tempera¬ 
tures. The dotted lines refer to the analogous results of Ref. [22], which also compute 
the rate for bremsstrahlung on nuclei. The effect of A Q to pair annihilation can be 
appreciated in this logarithmic scale only for the largest temperature 10 10 °K. 
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Fig. 6 . The regions in the T—p/p, e plane where each of the processes i)-iv) contribute 
for more than 90% to the total energy loss rate. We also show the contours for the rel¬ 
ative correction AQ/Q^ ot (see text) for the values 1%, 0%, —1%, —2%, —3%, —4%. 
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